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Abstract 

Using the effective Lagrangian method, we analyze the electroweak corrections to the anomalous 
dipole moments of lepton from some special two- loop diagrams where a closed neutralino/chargino 
loop is inserted into relevant two Higgs doublet one-loop diagrams in the minimal supersymmetric 
extension of the standard model with CP violation. Considering the translational invariance of 
loop momenta and the electromagnetic gauge invariance, we get all dimension 6 operators and 
derive their coefficients. After applying equations of motion to the external leptons, we obtain 
the anomalous dipole moments of lepton. The numerical results imply that there is parameter 
space where the contributions to the muon anomalous dipole moments from this sector may be 
significant. 
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I. INTRODUCTION 



At both aspects of experiment and theory, the magnetic dipole moments (MDMs) 
of lepton draw the great attention of physicists because of their obvious importance. The 
anomalous dipole moments of lepton not only can be used for testing loop effect in the 
standard model (SM), but also provide a potential window to detect new physics beyond 
the SM. The current experimental result of the muon MDM is 
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From the theoretical point of view, contributions to the muon MDM are generally divided 
into three sectors [2,3]: QED loops, hadronic contributions and electroweak corrections. The 
largest uncertainty of the SM prediction originates from the evaluation of hadronic vacuum 
polarization and light-by- light corrections. Depending on which evaluation of hadronic vac- 
uum polarization is chosen, the differences between the SM predictions and experimental 
result are given as 
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For the convenience of numerical discussion, we will adopt the second value in Eq|2j Within 
three standard error deviations, this difference implies that the present experimental data 
can tolerate new physics correction to the muon MDM as 

2.6 x 10~ 10 < Aa* p < 58.4 x 10" 10 . (3) 

In fact, the current experimental precision (6 x 10~ 10 ) already puts very restrictive bounds 
on new physics scenarios. In the SM, the electroweak one- and two-loop contributions 
amount to 19.5 x 10~ 10 and —4.4 x 10 -10 respectively. Comparing with the standard elec- 
troweak corrections, the electroweak corrections from new physics are generally suppressed 
by A^ w /A 2 , where A EW denotes the electroweak energy scale and A denotes the energy scale 
of new physics. 



Supersymmetry (SUSY) has been considered as a most prospective candidate for new 
physics beyond the SM. In the minimal supersymmetric extension of the SM (MSSM) with 
CP conservation, the supersymmetric one-loop contribution is approximately given by 

Aaf ~ 13 x lO- 10 (i^X) 2 tan/5sign(/iJ, (4) 

when all supersymmetric masses are assumed to equal a common mass A, and tan (3 = 
V2/V1 ^> 1. Where V\ and v 2 are the absolute values of the vacuum expectation values 
(VEVs) of the Higgs doublets and fi H denotes the //-parameter in the superpotential of 
MSSM. It is obvious that the supersymmetric effects can easily account for the deviation 
between the SM prediction and the experimental data. 

Actually, the two-loop electroweak corrections to the anomalous dipole moments of lepton 
are discussed extensively in literature. Utilizing the heavy mass expansion approximation 
(HME) together with the corresponding projection operator method, Ref. 4| has obtained 
the two-loop standard electroweak correction to the muon MDM which eliminates some 
of the large logarithms that were incorrectly kept in a previous calculation [||. Within 
the framework of MSSM with CP conservation, the authors of Ref. present the 

supersymmetric corrections from some special two-loop diagrams where a close chargino 
(neutralino) loop or a scalar fermion loop is inserted into those two-Higgs-doublet one-loop 
diagrams. Ref. 8j discusses the contributions to the muon MDM from the effective vertices 
H W^j, /io(-f^o)77 which are induced by the scalar quarks of the third generation in the 
CP conserving MSSM. 

In this paper, we investigate the electroweak corrections to the anomalous dipole moments 
of lepton from some special two- loop diagrams where a closed neutralino/chargino loop is 
inserted into relevant two Higgs doublet one-loop diagrams in the CP violating MSSM 
(FigHJ). Since the masses of those virtual fields {W ± , Z gauge bosons, neutral and charged 
Higgs, as well as neutralinos and charginos) are much heavier than the muon mass m^, we 
can apply the effective Lagrangian method to get the anomalous dipole moments of lepton. 
After integrating out the heavy freedoms mentioned above and then matching between the 
effective theory and the full theory, we derive the relevant higher dimension operators as 
well as the corresponding Wilson coefficients. The effective Lagrangian method has been 
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FIG. 1: Some two- loop self energy diagrams which lead to the lepton MDMs and EDMs in CP 
violating MSSM, the corresponding triangle diagrams are obtained by attaching a photon in all 
possible ways to the internal particles. In concrete calculation, the contributions from those mirror 
diagrams should be included also. 



adopted to calculate the two-loop supersymmetric corrections to the branching ratio of 
b — > S7 js|, neutron EDM 10] and lepton MDMs and EDMs 11]. In concrete calculation, we 
assume that all external leptons as well as photon are off-shell, then expand the amplitude of 
corresponding triangle diagrams according to the external momenta of leptons and photon. 
Using loop momentum translational invariance, we formulate the sum of amplitude from 
those triangle diagrams which correspond to the corresponding self-energy in the form which 
explicitly satisfies the Ward identity required by the QED gauge symmetry. Then we can get 
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all dimension 6 operators together with their coefficients. After the equations of motion are 
applied to external leptons, higher dimensional operators, such as dimension 8 operators, also 
contribute to the muon MDM and the electron EDM in principle. However, the contributions 
of dimension 8 operators contain an additional suppression factor mf/A 2 comparing with 
that of dimension 6 operators, where mi is the mass of lepton. Setting A ~ lOOGeV, one 
obtains easily that this suppression factor is about 10~ 6 for the muon lepton. Under current 
experimental precision, it implies that the contributions of all higher dimension operators 
(D > 8) can be neglected safely. 

We adopt the naive dimensional regularization with the anticommuting 7 5 scheme, where 
there is no distinction between the first 4 dimensions and the remaining D — 4 dimensions. 
Since the bare effective Lagrangian contains the ultraviolet divergence which is induced 
bv divergent subdiagrams, we give the renormalized results in the on-mass-shell scheme 
12j. Additional, we adopt the nonlinear gauge with £ = 1 for simplification 131 ] . This 
special gauge-fixing term guarantees explicit electromagnetic gauge invariance throughout 
the calculation, not just at the end because the choice of gauge-fixing term eliminates the 
, yW ± G Zf vertex in the Lagrangian. 

Within the framework of CP violating MSSM, the renormalization-group improved loop 
effects of soft CP violating Yukawa interactions related to scalar quarks of the third gen- 
eration cause the strong mixing among CP-even and CP-odd neutral Higgs. The linear 
expansions of the Higgs doublet H 1 and H 2 around the ground state are generally written 
as 
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where 9 is their relative phase. In the weak basis {(ft®, 
mass-squared matrix M 2 may be expressed as 
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V cos P 

Here, the concrete expressions of (A^|)y, {-M%p)ij 
the Higgs mass matrix Mjj is symmetric, we can diagonalize it by an orthogonal rotation 



sin P cos {3 

{M 2 SP ) 12 J 
can be found in the literature [14j . Since 



*Zj h dbS~. 



Z 1 M A Z H = diag(m 2 , m 2 , m 2 ). (7) 



H n 2 n 3 



Because of this strong mixing among the neutral Higgs, the couplings involving neutral Higgs 
are modified drastically comparing with that in CP conservating MSSM. Certainly, some 
diagrams in Figfj] have been discussed in Ref. [7\ where the authors apply the projecting 
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within the framework of CP 



operators to get the lepton MDMs (Eq.8~Eq.10 in Ref. 
conservating MSSM. On the other hand, the fermion electric dipole moments (EDMs) also 
offer a powerful probe for new physics beyond the Standard Model (SM). In the SM, the 
EDMs of leptons are fully induced by the CP phase of the Cabibbo-Kobayashi-Maskawa 



(CKM) matrix elements and they are predicted to be much smaller 15j than the present 



experimental precision Il7l | and beyond the reach of experiments in the near future. As 
for the MSSM, there are many new sources of the CP violation that can result in larger 
contributions to the EDMs of electron and neutron [18|, H9J. Taking the CP phases with 
a natural size of 0(1), and the supersymmetry mass spectra at the TeV range, we can 
find that the theoretical predictions on the electron and neutron EDMs at one-loop level 
already exceed the present experimental upper bound. In order to make the theoretical 
prediction consistent with the experimental data, one can generally adopt three approaches. 
One possibility is to make the CP phases sufficiently small, i.e. < 1CT 2 [18(]. One can 
also assume a mass suppression by making the supersymmetry spectra heavy, i.e. in the 
several TeV range [19], or invoke a cancellation among the different contributions to the 
fermion EDMs [20] . Since the lepton EDM is an interesting topic in both theoretical and 
experimental aspects, we as well present the lepton EDM by keeping all possible CP violating 
phases. 

This paper is composed by the sections as follows. In section [Til we introduce the effective 
Lagrangian method and our notations. Then we will demonstrate how to obtain the super- 
symmetric two-loop corrections to the lepton MDMs and EDMs. Section II III is devoted to 
the numerical analysis and discussion. In section IIVI we give our conclusion. Some tedious 
formulae are collected in appendix. 
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II. NOTATIONS AND TWO-LOOP SUPERSYMMETRIC CORRECTIONS 



The lepton MDMs and EDMs can actually be expressed as the operators 

C = 6 a Jct^I F 

£ BDM = -~ d t h^ 5 l . (8) 

Here, cr^ v — 7^/2, 1 denotes the lepton fermion, F is the electromagnetic field strength, 
m l is the lepton mass and e represents the electric charge, respectively. Note that the lepton 
here is on-shell. 

In fact, it is convenient to get the corrections from loop diagrams to lepton MDMs and 
EDMs in terms of the effective Lagrangian method, if the masses of internal lines are much 
heavier than the external lepton mass. Assuming external leptons as well as photon are all 
off-shell, we expand the amplitude of the corresponding triangle diagrams according to the 
external momenta of leptons and photon. After matching between the effective theory and 
the full theory, we can get all high dimension operators together with their coefficients. As 
discussed in the section HJ it is enough to retain only those dimension 6 operators in later 
calculations: 



1 



W) 

■ ' (4tt) 2 
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(9) 



with T> ^ = <9 M + ieA^ and u T = (1 7s)/2. When the equations of motion are applied to 
the incoming and outgoing leptons separately, only the operators actually contribute 
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to the MDMs and EDMs of leptons. We will only present the Wilson coefficients of the 
operators in the effective Lagrangian in our following narration because of the reason 
mentioned above. 

If the full theory is invariant under the combined transformation of charge conjugation, 
parity and time reversal (CPT), the induced effective theory preserves the symmetry after 
the heavy freedoms are integrated out. The fact implies the Wilson coefficients of the 
operators Of 3 satisfying the relations 

c* = cr, c+ = C e * , (10) 

where Cf (i = 1, 2, ■ ■ • , 6) represent the Wilson coefficients of the corresponding operators 
Of in the effective Lagrangian. After applying the equations of motion to the external 
leptons, we find that the concerned terms in the effective Lagrangian are transformed into 

c%o* + croj + c+o+ + cpo~ 

=► (C+ + C 2 * + Ct)Of + (C+* + C 2 + C?)Or 

= -^{^(C 2 + + C 2 * + Ct) 1<T I + i$(C+ + C 2 * + C+) To^is l}F, u . (11) 

Here, dt(- ■ •) denotes the operation to take the real part of a complex number, and 3(- • •) 
denotes the operation to take the imaginary part of a complex number. Applying Eq.® 
and Eq. lfTTT) . we finally get 

AQ t m 2 , , 
ttl = l^t^ + ° 2 * + C " ] ' 

di = - 2 -^Q(C+ + C 2 -* + C+). (12) 

In other words, the MDM of lepton is proportional to real part of the effective coupling 
C 2 + C 2 * + Cq, as well as the EDM of lepton is proportional to imaginary part of the 
effective coupling C 2 + C 2 * + C§ . 

Using the effective Lagrangian method, we present the one-loop supersymmetric contri- 



bution to muon MDM in ll| which coincides with the previous result in literature. Since 
the complication of analysis at two-loop order, we will adopt below a terminology where, for 
example, the "7/1^" contribution means the sum of amplitude from those triangle diagrams 
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(indeed three triangles bound together), in which a closed fermion (chargino/neutralino) 
loop is attached to the virtual Higgs and photon fields with a real photon attached in all 
possible ways to the internal lines. Because the sum of amplitude from those "triangle" di- 
agrams corresponding to each "self-energy" obviously respects the Ward identity requested 
by QED gauge symmetry, we can calculate the contributions of all the "self-energies" sep- 
arately Taking the same steps which we did in our earlier works 2I, lo|, llj], we obtain the 
effective Lagrangian that originates from the self energy diagrams in FigHJ In the bare ef- 
fective Lagrangian from the 'WW and 'ZZ' contributions, the ultraviolet divergence caused 
by divergent sub-diagrams can be subtracted safely in on-mass-shell scheme 12j. Now, we 
present the effective Lagrangian corresponding to the diagrams in Figfj] respectively. 



A. The effective Lagrangian from 7/1^ 



1,2,3) and 7G0 sector 



As a closed chargino loop is attached to the virtual neutral Higgs and photon fields, a 
real photon can be emitted from either the virtual lepton or the virtual charginos in the self 
energy diagram. When a real photon is emitted from the virtual charginos, the corresponding 
"triangle" diagrams belong to the typical two- loop Bar-Zee-type diagrams [21]. Within the 
framework of CP violating MSSM, the contributions from two-loop Bar-Zee-type diagrams 



221. When a 



to the EDMs of those light fermions are discussed extensively in literature 
real photon is attached to the internal standard fermion, the correction from corresponding 
triangle diagram to the effective Lagrangian is zero because of the Furry theorem, this point 
is also verified through a strict analysis. The corresponding effective Lagrangian from this 
sector is written as 



C 



2^2(4?r) 2 s 2 A 2 cos/3 



mi){-f) 1/2 ux h 
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with 



n% = (uluu L ) 2j (z H ) lk + (ul) i2 (u L )^z H ) 2k , 

•4 = ((C/t) ii (C/J 2j .sin/3+(C/t) i2 ( C /J ij cos/3)(Zj 3fc , = 1,2) . (14) 

Where the two unitary matrices U L R denote the left- and right-mixing matrices of charginos, 
A denotes the energy scale of new physics, and X{ = mf/A 2 respectively. We adopt the 
abbreviations: c = cos0, s w = sin0 , where is the Weinberg angle. The concrete 
expressions of Ti^ can be found in appendix. 

Accordingly, the lepton MDMs and EDMs from 7/1*. sector are written as 

hk _ V2e*Q f m%Z H ) lk ( x 4 y/2 

V2e 4 Q f m*(Z H ) 3k ttmf3 { x x f,i/2 

(AnY si A 2 V Xw J 2 ^ ' x x± ' X x f } ' 

d? fc * = - ^ Q f m '^ Q(Hti( X -£-) 1/2 T 2 (x h ,x ± ,x J 
1 v / 2(47r) 4 s 2 A 2 cos /3 1 *^ w fc 

which are enhanced by large tan/5. Note here that the corrections from this sector to the 
MDM of lepton depend on a linear combination of real parts of the effective couplings 
and A%, and the corrections from this sector to the EDM of lepton depend on a linear 
combination of imaginary parts of the effective couplings TL^ and A\. In the limit x ± x hk , 
the above expressions can be simplified as 

ai (47r) 4 S 2 A 2 cos/3 M[nj y x > * ± ™ ± dx 4 

_ ^g / m 2 (Zj 3fc tan/j f ^y 2[ ln^ 9 

(4vr) 4 S 2 A 2 ^ AJ l x J i x . + - i™ ± , W 



± ' 3 



J**- ^Q f mXZ H ) lk ^k/xj^Mx^ dj x x ) 



v / 2(4tt) 4 s 2 A 2 cos /3 " v a; w 1 1 x^ ± x x ±^ x x ± dx^ 

• 5 Q fmi (Z H ) 3k tan/3 fc A 

V2(4vr) 4 s 2 A 2 ^ l iJ ^ w y *J^> x ±dx 



e»q^(Z ) tt tan/3 ^ )Va 
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Similarly, we can formulate the corrections from 7G0 sector to the effective Lagrangian 



as 



with 



x 



ClG 2V2(4tt) 2 s 2 A 2 



K(^)(^) 1/2 T 2 (x, 
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iZ(B u ) v \ (x z , x , x ) (0+ - O;) \ , 

X w * * ) 



(17) 



By = -W) a (UX cos (3 + (Ul) a (UX sM , = 1, 2) 



(18) 



Correspondingly, the corrections to the lepton MDMs and EDMs from this sector are: 



lG V2e 4 Q f m* 
(4tt) 4 s 2 A 2 
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^B ti )(^) 1/2 T 2 (x z ,x ± , 
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^ 7 g _ e 5 Q f m l 
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9f(B«)(— J Ti(x z ,x ± ,x J . 



(19) 



v / 2(4tt) 4 s 2 A 2 

The corrections from this sector to the MDM of lepton are proportional to real parts of 
the effective couplings Bu, and the corrections from this sector to the EDM of lepton are 
proportional to imaginary parts of the effective couplings Bu, separately. In the limit x , 
x z , we have 



y/2e 4 Q f m* 
(4tt) 4 s 2 A 2 



M(Bii)C—) 
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xf x± "*„± 



v / 2(4tt) 4 s 2 A 2 



) lim 



dx 



x ± ^x ± ^ ± 

X J X i X. 



-(^i(x ± ,x ± ) . 



(20) 



Using the concrete expression of <pi(x,y) presented in appendix, one can verify easily 
that the corrections to the lepton MDMs and EDMs from the sectors are suppressed by the 
masses of charginos asm ± > m h , m z (i — 1, 2). 

X,- k 



B. The effective Lagrangian from Zhk (ZGq) sector 



As a closed chargino loop is attached to the virtual Higgs and Z gauge boson fields, a 
real photon can be attached to either the virtual lepton or the virtual charginos in the self 
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energy diagram. When a real photon is attached to the virtual lepton, the corresponding 



amplitude only modifies the Wilson coefficients of the operators Of in the effective La- 
grangian after the heavy freedoms are integrated out. In other words, this triangle diagram 
does not contribute to the lepton MDMs and EDMs. A real photon can be only attached to 
the virtual lepton as the closed loop is composed of neutralinos, the corresponding triangle 
diagram does not affect the theoretical predictions on the lepton MDMs and EDMs for the 
same reason. Considering the points above, we formulate the contributions from Zh$ sector 
to the effective Lagrangian as 
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16V2(47r) 2 s^c2Q / A2 cos/3 
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(21) 



with 



eS = 2^cos20 w + (C/t) ii ([/J i . , 

€5 = 2fycos20 w + (17 K ) M , = 1,2) , (22) 

where the concrete expressions of the functions ^ .(rri,^), ^1,2(^1, £2, ^3, Xi) are listed in 
appendix. Additional, T z is the isospin of lepton, and Q f is the electric charge of lepton, 
respectively. Using Eqf2Tl we get the corrections to the lepton MDMs and EDMs from Zhk 
sector as 



Zh k 
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(23) 



The above equations contain the suppression factor 1— 4s^ because Q f — — 1 and = — 1/2 
for charged leptons. The corrections from this sector to the MDM of lepton are decided by a 
linear combination of real parts of the effective couplings 7i k £ L + TC ^£ R and A k £ —A ' E 1 *, 
and the corrections from this sector to the EDM of lepton are decided by a linear combination 
of imaginary parts of the effective couplings 7i k £ L — Ti !<£ R and A k £ L + A k, ^E R . In the limit 
x ± , x ± x z , x h , Eqj23]can be approximated as 
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Similarly, the contribution from ZGq sector to the effective Lagrangian is 

4 r x 
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and the contributions to the lepton MDMs and EDMs are: 
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(26) 



Here, the corrections from this sector to the MDM of lepton are proportional to real parts 
of the effective couplings B £ L — B^ £ R , and the corrections from this sector to the EDM of 
lepton are proportional to imaginary parts of the effective couplings B-.£ L . + B\.£^. When 
x ± , x ± ^> x z , Eql26] can be approached by 
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C. The effective Lagrangian from 7Z sector 



When a closed chargino loop is attached to the virtual 7 and Z gauge bosons, the 
corresponding correction to the effective Lagrangian is very tedious. If we ignore the terms 
which are proportional to the suppression factor 1 — 4s 2 , the correction from this sector to 
the effective Lagrangian is drastically simplified as 



C 



-,z 



UL_tR\ lim T 

V*" / x ^x , " 
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X,- x . 



8(4tt) 2 s 2 c 2 A 2 
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(28) 

Using the definitions of the matrices in Eq. (l22p . one can find that the effective couplings 
£ — £ (i = 1,2) are real. Correspondingly, the correction to the lepton MDMs from this 
sector is written as 

,2 
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17T l 4 .S^ H A 2 \ " " / X , ^X , XT XT 



4(4vr) 4 s 2 c 2 A 2 

V / WW 



x i x i 



(29) 



and the correction to the lepton EDMs is zero. In the limit x ± ^> x z , we can approximate 



the correction to the lepton MDMs from this sector as 



z e 4 Q f m 2 



4(4tt) 4 s 2 c 2 A 2 



(e - e 
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10 In x , — 2 In x 

lO x 

+ - 



3x 



+ lim [2x 



x i x i 



t dx 2 , 



_j_ ^ ^-j 



(30) 



D. The effective Lagrangian from W^H^ (W^G^) sector 



As a closed chargino-neutralino loop is attached to the virtual W gauge boson and 
charged Higgs H T , the induced Lagrangian can be written as 

4 j n ( x 



e tan (3 



16(47r)Vc w ^A 2 l^w ; 31 
- cos pgf(«) O- + ( sin f3&% ((% - cos W% ((%)0: 



j /■>•* ry ry rr> 



+ - 

v x„ 



v 1/2 



) Tf sin ^ 1L c fi - cos pg 1R ( L ) o: 
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+ ( sin (3(g 1L ) t ((% - cos P(G 1R )l (C%)o: 



x 



+ (^) F5(x w ,^ ± ,x xO ,X x± )[(sin/3^C; + cos/3^C;)0 

w * 3 

+ (^) V2 ^ w , x ff± , x^ , x^ ± ) [( sin + cos 1lJ]\'-) 0' 6 

+ ( sin ^% ((% + cos (,'■) )cr 



(31) 



with 



cs=^j(^)„-^i(^j !y , 

C? = ^(^), + ^AT 3 ^t) j2 , 
(* = !,•••, 4, j = 1,2) . 



(32) 



Here, the 4x4 matrix A/" denotes the mixing matrix of the four neutralinos (« = 1, • • • , 4) . 
The corresponding corrections to the lepton MDMs and EDMs are respectively expressed 

as 



aY H = 



{ {-£) 1/2 Fs(x w , x H± , x x , , x xf )»( sin >'(;"(/ - cos /^C, 



e 4 m 2 tan /? 
■4(4vr)%4 CwA 2 

+ (-f) V V 4 (x w , x„ ± , x x0 , x )»( sin iG"x> ! - cos /^Cj) 



x 



X 
X 



)9fc( sin^ 1L C L + cos/^C*) 



x? \ 1/2 



X 



F 6 (x v ,x H± ,X x0 ,X x± n(smpgl L C R + cospg^X) , 



dY H = 



e 5 m tan/3 



x 



8(4vr) 4 s 4 c w A 2 1 V x 



1/2 



F 3 (x w , x h± , x x0 , x ) Of ( sin ic;; \ '; - cos #"<<,") 



jSfsin^C^-cos^ 1 ^) 



+ (-^) 1/2 ^5 (x w , x ff± , x , x )3 ( sin + cos /3^CJ) 
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+ CfY /2 Fe(x w ,x H± ,x x0 ,x )»(sin/3^ + cos/3^C;) [ , (33) 



X... ' X i X J 



where the concrete expressions of ^3,4,5,6 can be found in appendix. The corrections from this 
sector to the MDM of lepton are decided by a linear combination of real parts of the effective 
couplings smf3Q 1L ( L - cos 8G 1R ( R , sin f3G 1L ( R - cos 3G 1R ( L , sin 3G 1L ( L + cos 3G 1R ( R } as 

A 3% 13 31 13 31 13 31 13 31 13 31 13 

well as sin 3Q 1L C R + cos 3Q XR C L ', and the corrections from this sector to the EDM of lepton 

31 13 31 %3 

are decided by a linear combination of imaginary parts of those effective couplings. Using 
the asymptotic expansion of the two-loop vacuum integral y, z) presented in appendix, 
we can simplify the expressions of Eqj33] in the limit x „ , x ± 3> x . 

As a closed chargino-neutralino loop is attached to the virtual W ± gauge boson and 
charged Goldstone G T , the corresponding corrections to the lepton MDMs and EDMs are 
similarly formulated as 



' ("f ) 11 V 3 (s w , x wt * x? , x xf )»( cos 8g]X + sin 



4( 



4(4tt) 4 s 4 / c w A 2 x 
5 x? \ 1/2 



+ (^) F 4 (x w ,x w ,X xr X x JK(cos^C; + sin/?^C; 
+ ("^) 1/2 ^5 (* w , x w , x x0 , x )»( cos 3QIX - sin 



X 
X 



+ (-r) 1/2 ^ . ^ , ^ x o , *J»( cos - sin *;'.v 



= - - e m > — J (J^) 1/2 F 3 (a; w> a; wJ a; ,x ± )9f( cos/3£ 1L C L + sin^V) 

V / w W \ w J 



+ {-ff' 2 F^x^x^x^{ cos fl^C? + sin 

w 1 3 

+ (■£■) 1/2 ^ (x w , x w , ^ , x x , ) 9f ( cos kV'x 1 / - sin /?SfC*) 

w 1 3 ) 

Similarly, the corrections from this sector to the MDM of lepton depend on a linear combina- 
tion of real parts of the effective couplings cos 3Q 1L ( L + sm 3Q 1R ( R , cos 3Q 1L ( R +sm 3Q lR C L , 
cos 3Q lL C L — sin 3Q 1R C R , as well as cos 3Q lL C R — sin BQ lR C L , and the corrections from this 

ji ij 3% ij ji ij ji ij 

sector to the EDM of lepton depend on a linear combination of imaginary parts of those 
effective couplings. 
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The contributions from those above sectors to effective Lagrangian do not contain ul- 
traviolet divergence. In the pieces discussed below, the coefficients of high dimensional 
operators in effective Lagrangian contain ultraviolet divergence that is caused by the di- 
vergent subdiagrams. In order to obtain physical predictions of lepton MDMs and EDMs, 
it is necessary to adopt a concrete renormalization scheme removing the ultraviolet diver- 
gence. In literature, the on-shell renormalization scheme is adopted frequently to subtract 
the ultraviolet divergence which appears in the radiative electroweak corrections [12]. As 
an over-subtract scheme, the counter terms include some finite terms which originate from 
those renormalization conditions in the on-shell scheme beside the ultraviolet divergence 
to cancel the corresponding ultraviolet divergence contained by the bare Lagrangian. In 
the concrete calculation performed here, we apply this scheme to subtract the ultraviolet 
divergence caused by the divergent subdiagrams. 



E. The effective Lagrangian from the ZZ sector 

The self energy of Z gauge boson composed of a closed chargino loop induces the 
ultraviolet divergence in the Wilson coefficients of effective Lagrangian. Generally, the un- 
renormalized self energy of the weak gauge boson Z can be written as 



£i (p) = A 2 A*g, u + [A* + j- 2 Al ) {p 2 g, u - p,p v ) + (B{ + ^flf ) PflPv . (35) 
Correspondingly, the counter terms are given as 

K?(p) = -(<H 2 + m 1 SZ *)9^ - K(PV - PmP") • (36) 
The renormalized self energy is given by 

^ Z » = S z » + S^(p). (37) 

For on-shell external gauge boson Z, we have [l 



2 2 = ' 



KSpV(p) 

J% P ^ t > )eU{p) = e ^ h (38) 
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FIG. 2: The counter term diagram to cancel the ultraviolet caused by the self energy of Z boson. 



where e{p) is the polarization vector of Z gauge boson. From Eq. (138]) . we get the counter 
terms 

,2 



m 



5Z z =Al + J ^A z 2 = At + xA z 2 , 
Sm 2 = A*A 2 -m 2 5Z . 

z U z z 



(39) 



Accordingly, the effective Lagrangian originating from the counter term diagram (FigfSJ) 
can be formulated as 
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Here, e = 2 — D/2 with D representing the time-space dimension, and x H = A 2 E /A (A RE 
denotes the renormalization scale). 

As a result of the preparation mentioned above, we can add the contributions from the 
counter term diagram to cancel the corresponding ultraviolet divergence contained by the 
bare effective Lagrangian. Using the definitions of the matrices £ L,i? in Eq. fl22l . we derive 
tL _ cl* tR _ tR* ^he resu ited theoretical predictions on the lepton MDMs and EDMs 

ij ji ij ji 

are respectively written as 

9 /■ 
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{Tf-Qr 



Qy 



T s (x„,x 



(41) 



In other words, the corrections from this sector to the MDM of lepton are decided by a 
linear combination of the real effective couplings |£^| 2 ± l£^l 2 an< ^ rea ^ P ar ts of the effective 
couplings £, L £, R , and the corrections from this sector to the EDM of lepton are proportional 
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to imaginary parts of the effective couplings £, L !; R - 

l 3 3 l 

Because a real photon can not be attached to the internal closed neutralino loop, the 
corresponding effective Lagrangian only contains the corrections to the lepton MDMs: 



O, )X 



e A Q f mf 
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(42) 



A). 



(43) 



In order to get Eq. 042p . we apply unitary property of the matrices rj L,R . The corrections 
from this sector to the MDM of lepton depend on a linear combination of the real effective 
couplings \f]^\ 2 ± |?7^| 2 and real parts of the effective couplings vf-V^ an< ^ *he corrections 
from this sector to the EDM of lepton are proportional to imaginary parts of the effective 
couplings T] L j] R . 

^3 3 l 

We can also simplify Eq. (J4T]1 and Eq. (f4"2l in the limit x ± ,x ± ,x ,x 3>x z using the 

X i Xj X t Xj 

asymptotic expansion of $(x,y,z). The concrete expressions of T 4 ~ T 9 can be found in 
appendix. 



F. The effective Lagrangian from the WW sector 

Similarly, the self energy of W gauge boson composed of a closed chargino-neutralino 
loop induces the ultraviolet divergence in the Wilson coefficients of effective Lagrangian. 
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FIG. 3: The counter term diagram to cancel the ultraviolet caused by the self energy of W boson 
and electroweak radiative corrections to ^yW + W~ vertex. 



Accordingly, the unrenormalized W self energy is expressed as 



p- 



The corresponding counter terms are given as 

The renormalized self energy is given by 

£ w (p) = £ w (p) + £ wc (p) 
For on-shell external gauge boson W , we have 

£ w (p)e"(p) , , = 
lim 
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where e(p) is the polarization vector of W gauge boson. Inserting Eq. ( 1441) and Eq. (145p 
into Eq. ( jUj) , we derive the counter terms for the self energy as 

2 

772 

5m 2 = A™A 2 -m 2 5Z . (48) 

Differing from the analysis in the ZZ sector, we should derive the counter term for the vertex 
'yW + W~ here since the corresponding coupling is not zero at tree level. In the nonlinear 
gauge with £ = 1, the counter term for the vertex / yW + W~ is 



iSC lW + w - = ie ■ SZ^ 



g^u{h - k 2 ) p + g vp {k 2 - h)n + 9pit[fa- k i)y > ( 49 ) 



where ki (i = 1, 2, 3) denote the injection momenta of W ± and photon, and //, u, p denote 
the corresponding Lorentz indices respectively. 

We present the counter term diagrams to cancel the ultraviolet divergence contained in 
the bare effective Lagrangian from WW sector in Figj3l and we can verify that the sum of 
corresponding amplitude satisfies the Ward identity required by the QED gauge invariance 
obviously. Accordingly, the effective Lagrangian originating from the counter term diagrams 
can be written as 
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Finally, we get the renormalized effective Lagrangian from the WW sector: 
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(50) 
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Correspondingly, the resulted lepton MDMs and EDMs are respectively formulated as 
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(52) 



In a similar way, the corrections from this sector to the MDM of lepton depend on a linear 
combination of the real effective couplings \(^\ 2 ± |0j| 2 an d real parts of the effective cou- 
plings ( L ( R *, and the corrections from this sector to the EDM of lepton are proportional to 
imaginary parts of the effective couplings ( L ( R *. 



III. NUMERICAL RESULTS AND DISCUSSION 



With the theoretical formulae derived in previous section, we numerically analyze 
the dependence of the muon MDM and EDM on the supersymmetric parameters in the 
CP-violating scenario here. In particular, we will present the dependence of the muon 
MDM and EDM on the supersymmetric CP phases in some detail. In order to make the 
theoretical predictions on the electron and neutron EDMs satisfying the present experimental 
constraints, we adopt the cancelation mechanism among the different contributions to the 
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fermion EDMs 20] .Within three standard error deviations, the present experimental data 
can tolerate new physics correction to the muon MDM as 2.6 x 1CT 10 < Aa /t < 58.4 x 
1CT 10 . Since the neutralinos Xi (* = 1j 2, 3, 4) and charginos xt (* = 1> 2) appear as 
the internal intermediate particles in the two-loop diagrams which are investigated in this 
work, the corrections of these diagrams will be suppressed strongly when the masses of 
neutralinos and charginos are much higher than the electroweak scaleQ]. To investigate if 
those diagrams can result in concrete corrections to the muon MDM and EDM, we choose 
a suitable supersymmetric parameter region where the masses of neutralinos and charginos 
are lying in the range M x < 600 GeV. 

The MSSM Lagrangian contains several sources for CP violating phases: the phases of 
the /i-parameter in the superpotential and the corresponding bilinear coupling of the soft 
breaking terms, three phases of the gaugino masses, and the phases of the scalar fermion 
Yukawa couplings in the soft Lagrangian. As we do not consider the spontaneous CP 
violation in this work, the CP phase of soft bilinear coupling vanishes due to the neutral 
Higgs tadpole conditions. Additional, the CP violation would cause changes to the neutral- 
Higgs-quark coupling, the neutral Higgs-gauge-boson coupling and the self-coupling of Higgs 
boson. A direct result of above facts is that no absolute limits can be set for the Higgs bosons 
masses from the present combined LEP data 23J]. For security, we take the lower bound 
on the mass of the lightest Higgs boson as m h > 60GeV [14j in the numerical analysis. In 
order to obtain the mixing matrix of neutral Higgs in CP violating MSSM, we include the 
subroutine fillhiggs.f from the Package CPsuperH [24| in our numerical code. Furthermore, 
we take the pole mass of top quark m t (pole) = 175 GeV, the pole mass of charged Higgs 
m .(pole) = 300 GeV, the running masses m&(mt) = 3 GeV, m T {m t ) = 1.77GeV, the mass 
parameters of scalar fermions in soft terms as m - = m~ =m. =m. = m. = 500 GeV, 

r U 3 D a E 3 Q 3 L 3 

the Yukawa couplings of scalar fermions as \A t \ = \Ab\ = \A T \ = 1 TeV and 4>A t = 4>A b = 
(f>A T = 7r/2. Fixing above parameters and assuming tan/? > 3, we find that the mass of 
the lightest neutral Higgs is well above 115 GeV by scanning the parameter space of CP 
violating MSSM. In other words, one no longer worries about the constraint from Higgs 
sector with the above assumptions on the parameter space of CP violating MSSM. With no 
loss of generality, we also take the supersymmetric parameters \mi\ = \rri2\ = 500 GeV and 
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FIG. 4: The supersymmetric corrections to the muon MDM and EDM d„ vary with the CP 
violating phase (f) mi when = 200 GeV, <fi m2 = = and tan/3 = 10, 50, where the solid 
lines stand for the one-loop corrections with tan /3 = 10, the dot lines stand for the results in- 
cluding two-loop supersymmetric corrections with tan = 10; the dash lines stand for the one-loop 
corrections with tan/3 = 50, the dash-dot lines stand for the results including two-loop supersym- 
metric corrections with tan = 50. The gray band in diagram (a) is the region allowed by the g — 2 
experimental data within 3 standard errors. 



m. = m~ = A a /2 = 500 GeV in this work. 

Taking \fi H \ = 200 GeV, ma = M =0 and tan/3 = 10, 50, we plot the muon MDM a M 
and EDM versus the CP phase (f) m in FigJU As tan /3 = 10, the one-loop supersymmetric 
correction to the muon MDM (solid-line in FigJU^a)) reaches 7 x 10 -10 and can account 
for the deviation between the SM prediction and experimental data. Comparing with one- 
loop supersymmetric contribution, two-loop contribution depends on the supersymmetric 
parameters in a different manner. Including the two-loop corrections, the supersymmetric 
contribution to the muon MDM is modified about 10%. Since the gaugino mass mi 
affects the theoretical prediction only through the mixing matrix of neutralinos, the muon 
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MDM a M varies with the CP phase <p m (solid line for one-loop result and dot line for the 
result including two-loop corrections in FigH^a)) very mildly. Meanwhile the supersymmetric 
contribution to the muon EDM including two-loop corrections at the largest CP violation 
<f> mi = 7r/2 is still below 10 _24 e • cm (dot line FigJH(b)), and it is very difficult to observe the 
muon EDM of this level in next generation experiments with precision 10™ 24 e ■ cm [17]. As 
tan/3 = 50, one-loop supersymmetric correction to the muon MDM exceeds 35 x 1CT 10 
(dash line in FigJU(a)), and can ameliorate easily the discrepancy between the SM prediction 
and experiment. Because the dominant two-loop supersymmetric corrections originating 
from the jhk, W ± H Zf sectors are enhanced by large tan/3, the relative modification from 
two-loop supersymmetric corrections to one-loop result is 15% roughly (dash-dot line in 
FigHfa)). As for the muon EDM d^, one-loop supersymmetric result together with two-loop 
supersymmetric corrections are all enhanced by large tan/3. The contribution including 
two-loop supersymmetric corrections is well above 10~ 24 e ■ cm at the largest CP violation 
(f) m = tt/2, and it is hopeful to detect the muon EDM of this level in the near future. 

Taking \fi H \ = 200 GeV, <f> mi =<P liH =0 and tan/3 = 10, 50, we plot the muon MDM 
and EDM d^ versus the CP phase ma in FigJSJ As tan (3 = 10, the one-loop supersymmetric 
correction to the muon MDM (solid-line in FigJS^a)) always lies in the range ja^j < 8 x 10~ 10 
varying with the CP phase </> ma . The relative modification from the two-loop supersymmetric 
corrections to the one-loop prediction is below 5% when tan/3 = 10. Since the gaugino 
mass m 2 affects the theoretical prediction through the mixing matrices of neutralinos and 
charginos simultaneously, the muon MDM a M depends on the CP phase m2 (solid line 
for one-loop result and dot line for the result including two- loop corrections in Figj5]^a)) 
strongly. Meanwhile the supersymmetric contribution to the muon EDM including two-loop 
corrections at the largest CP violation </> ma = 7r/2 is about 10 _23 e • cm (dot line Figj5 b)) 
which can be observed in next generation experiments with precision 10~ 24 e ■ cm 17| . 
When tan/3 = 50, one-loop supersymmetric correction to the muon MDM a M is enhanced 
drastically. Because the dominant two-loop supersymmetric corrections originating from 
the 7/1*;, W ± H T sectors are also enhanced by large tan/3, the relative modification from 
two-loop supersymmetric corrections to one-loop result is 15% roughly (dash-dot line in 
FigJS^a)). As for the muon EDM d^, one-loop supersymmetric result together with two-loop 
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FIG. 5: The supersymmetric corrections to the muon MDM and EDM d„ vary with the CP 
violating phase 4> m2 when = 200 GeV, 4> mi = = an d tan/3 = 10, 50, where the solid 
lines stand for the one-loop corrections with tan (3 = 10, the dot lines stand for the results in- 
cluding two-loop supersymmetric corrections with tan = 10; the dash lines stand for the one-loop 
corrections with tan/3 = 50, the dash-dot lines stand for the results including two-loop supersym- 
metric corrections with tan = 50. The gray band in diagram (a) is the region allowed by the g — 2 
experimental data within 3 standard errors. 



supersymmetric corrections are all enhanced by large tan (3. The contribution including two- 
loop supersymmetric corrections at the largest CP violation ma = 7r/2 is about 4x 10~ 23 e-cm 
which can be detected easily in next generation experiments. 

Taking \n H \ = 200 GeV, <fi = <f) m2 = and tan/3 = 10, 50, we plot the muon 
MDM a M and EDM versus the CP phase in Figj6j As tan/3 = 10, the one-loop 
supersymmetric correction to the muon MDM (solid-line in Figj6](a)) always lies in the range 
< 8 x 10~ 10 varying with the CP phase • The relative modification from the two- 
loop supersymmetric corrections to the one-loop prediction is below 5% when tan/3 = 10. 
Since the /i parameter fi H affects the theoretical prediction through the mixing matrices 
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FIG. 6: The supersymmetric corrections to the muon MDM and EDM d„ vary with the CP 
violating phase </> when |/x H | = 200 GeV, <j> mi = <^ ma = and tan/3 = 10, 50, where the solid 
lines stand for the one-loop corrections with tan (3 = 10, the dot lines stand for the results in- 
cluding two-loop supersymmetric corrections with tan = 10; the dash lines stand for the one-loop 
corrections with tan/3 = 50, the dash-dot lines stand for the results including two-loop supersym- 
metric corrections with tan = 50. The gray band in diagram (a) is the region allowed by the g — 2 
experimental data within 3 standard errors. 



of neutralinos and charginos simultaneously, the muon MDM a M varies with the CP phase 
M (solid line for one-loop result and dot line for the result including two-loop corrections 
in Figj6](a)) drastically. Meanwhile the supersymmetric contribution to the muon EDM 
including two-loop corrections at the largest CP violation </> = 7r/2 is below 10 _23 e • cm 
(dot line FigJB^b)). Because the dominant two-loop supersymmetric corrections originating 
from the jhk, W ± H T sectors are enhanced by large tan/3, the relative modification from 
two-loop supersymmetric corrections to one-loop result is 15% roughly (dash-dot line in 
FigJB^a)) at CP conservation when tan/3 = 50. One-loop supersymmetric correction to the 
muon EDM is enhanced by large tan j3. Comparing with one-loop contribution, two-loop 
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FIG. 7: The supersymmetric corrections to the muon MDM vary with the //-parameter /t H 
when (f> mi = 4> m2 = (j)^ = and tan j3 = 20, 50, where the solid lines stand for the one-loop 
corrections with tan (5 = 20, the dot lines stand for the results including two-loop supersymmetric 
corrections with tan = 20; the dash lines stand for the one-loop corrections with tan/3 = 50, the 
dash-dot lines stand for the results including two-loop supersymmetric corrections with tan = 50. 
The gray band is the region allowed by the g — 2 experimental data within 3 standard errors. 



corrections are negligible. The contribution including two-loop supersymmetric corrections 
is about 4 x 10~ 23 e • cm, which can be detected in the near future 17]. 



Taking tan/5 = 20, 50 and <p mi = </> m2 = </> M ^ = 0, we plot the muon MDM a M versus the 
//-parameter fi H in FigjTl The gray band is the region allowed by present experimental data 
within 3 standard errors. Because the supersymmetric corrections to the muon MDM a M are 
negative for n H < 0, the corresponding parameter space is already ruled out by the present 
g — 2 experimental data. Comparing with the one-loop supersymmetric results (solid line for 
tan/3 = 20 and dash line for tan/3 = 50 respectively), the contributions including two- loop 
supersymmetric corrections are enhanced about 15% when /i H = 150 GeV. Along with the 
increasing of n H , the two-loop corrections become more and more trivial. 
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IV. CONCLUSIONS 



In this work, we analyzed the two-loop supersymmetric corrections to the muon MDM 
and EDM by the effective Lagrangian method in the CP violating MSSM. In the concrete 
calculation, we keep all dimension 6 operators. The ultraviolet divergence caused by the 
divergent sub-diagrams is removed in the on-shell renormalization schemes. After apply- 
ing the equations of motion to the external leptons, we derive the muon MDM and EDM. 
Numerically, we analyze the dependence of the muon MDM and EDM on supersym- 
metric CP violating phases. As discussed above, is decided by real parts of the effective 
couplings, and is decided by imaginary parts of the effective couplings after the heavy 
freedoms are integrated out. Adopting our assumptions on parameter space of the MSSM 
and choosing tan/5 = 50, we find that the correction from those two-loop diagrams to a M is 
4 x 10~ 10 roughly for the case of CP conservation, which lies in the order of present exper- 
imental precision in magnitude. In other words, the present experimental data put a very 
restrictive bound on the real parts of those effective couplings. Additional, the contribution 
to d^ from this sector is sizable enough to be experimentally detected with the experimental 
precision of near future. 
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APPENDIX A: THE FUNCTIONS 

We list the tedious expressions of the functions adopted in the text 

x l In- 7 x — y l In- 7 y 

,3 x — y 
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The concrete expression of $(x, y, z) can be found in lCj, |25|. In the limit z <C x,y, we 
can expand &(x,y,z) according z as 
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Here, the function Q(x,y) is defined as 



2 

Q(x,y) =lnxln^-21n(x-y)ln^-2L« 2 (^) + — . (A8) 
x x x 3 
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